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I. INTRODUCTION 



The issue to extend General Relativity (GR) to alternative theories of gravity has recently become dramatically 
urgent due to the missing matter problem at all astrophysical scales and the accelerating behavior of cosmic fluid, 
detected by SupcrNovae la used as standard candles. Up to now, no final answer on new particles has been given at 
fundamental level so Dark Energy and Dark Matter constitute a puzzle to be solved in order to achieve a self-consistent 
picture of the observed Universe, /(i?)- gravity, where f(R) is a generic function of the Ricci scalar R, comes into the 
game as a straightforward extension of GR where further geometrical degrees of freedom are considered instead of 
searching for new material ingredients From an epistemological point of view, the action of gravity is not selected 
£h . a priori, but it could be "reconstructed", in principle, by matching consistently the observations This approach 

can be adopted considering any function of the curvature invariants as R^ V R^ V , ROR and so on. 

From a genuine mathematical point of view, alternative theories of gravity pose the problem to recover or extend 
CN ' the well-established results of GR as the initial value problem 0, the stability of solutions and, in particular, the issue 
^ of finding out new solutions. As it is well known, beside cosmological solutions, spherically and axially symmetric 
■ solutions play a fundamental role in several astrophysical problems ranging from black holes to active galactic nuclei. 

Alternative gravities, to be consistent with results of GR, should comprise solutions like Schwarzschild and Kerr ones 
. but present, in general, new solutions that could be physically interesting. Due to this reason, methods to find out 
exact and approximate solutions are particularly relevant in order to check if observations can be framed in Extended 
Theories of Gravity . 

Recently, the interest in spherically symmetric solutions of /(i?)-gravity is growing up. In 0, solutions in vacuum 
\ have been found considering relations among functions that define the spherical metric or imposing a costant Ricci 
curvatue scalar. The authors have reconstructed the form of some /(i?)-models, discussing their physical relevance. In 
Q, the same authors have discussed static spherically symmetric solutions, in presence of perfect fluid matter, adopting 
the metric formalism. They have shown that a given matter distribution is not capable of globally determining the 
$H ' functional form of f(R). Others authors have discussed in details the spherical symmetry of /(i?)-gravity considering 
also the relations with the weak field limit. Exact solutions are obtained for constant Ricci curvature scalar and for 
Ricci scalar depending on the radial coordinate. In particular, it can be considered how to obtain results consistent 
with GR assuming the well-known post-Newtonian and post-Minkowskian limits as consistency checks. @. 

In this paper, we want to seek for a general metod to find out axially symmetric solutions by performing a complex 
coordinate transformation on the spherical metrics. Since the discovery of the Kerr solution (loj . many attempts 
have been made to find a physically reasonable interior matter distribution that may be considered as its source. 
For a review on these approaches see [ill, IH|- Though much progress has been made, results have been generally 
disappointing. As far as we know, nobody has obtained a physically satisfactory interior solution. This seems 
surprising given the success of matching internal spherically symmetric solutions to the Schwarzschild metric. The 
problem is that the loss of a degree of symmetry makes the derivation of analytic results much more difficult. Severe 
restrictions are placed on the interior metric by maintaining that it must be joined smoothly to the external axially 
symmetric metric. Further restrictions are placed on the interior solutions to ensure that they correspond to physical 
objects. 

Furthermore since the axially symmetric metric has no radiation field associated with it, its source should be also 
non-radiating. This places even further constraints on the structure of the interior solution [l3j . Given the strenuous 
nature of these limiting conditions, it is not surprising to learn that no satisfactory solution to the problem of finding 
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sources for the Kerr metric has been obtained. In general, the failure is due to internal structures whose physical 
properties are unknown. This shortcoming makes hard to find consistent boundary conditions. 

Newman and Janis showed that it is possible to obtain an axially symmetric solution (like the Kerr metric) by 
making an elementary complex transformation on the Schwarzschild solution [l5j |. This same method has been used 
to obtain a new stationary and axially symmetric solution known as the Kerr-Newman metric [l6j . The Kerr- Newman 
space-time is associated to the exterior geometry of a rotating massive and charged black-hole. For a review on the 



Newman- Janis method to obtain both the Kerr and Kerr-Newman metrics see 17 1. 

By means of very elegant mathematical arguments. Schiffcr et al. fl8| have given a rigorous proof to show how the 
Kerr metric can be derived starting from a complex transformation on the Schwarzschild solution. We will not go 
into the details of this demonstration, but point out that the proof relies on two main assumptions. The first is that 
the metric belongs to the same algebraic class of the Kerr-Newman solution, namely the Kerr-Schild class [19 ]. The 
second assumption is that the metric corresponds to an empty solution of the Einstein field equations. In the case we 
are going to study, these assumptions are not considered and hence the proof in pj} is not applicable. It is clear, by 
the generation of the Kerr-Newman metric, that all the components of the stress-energy tensor need to be non-zero 
for the Newman- Janis method to be successful. In fact, Giirses and Giirsey, in 1975 [20(, showed that if a metric can 
be written in the Kerr-Schild form, then a complex transformation "is allowed in General Relativity." In this paper, 
we will show that such a transformation can be extended to /(i?)-gravity. 

The outline of this paper is as follows. In the SecflTl we introduce the /(-R)-gravity action, the field equations and 
ive some general remarks on spherical symmetry. In Sec IIII1 we give a summary of the Noether Simmetry Approach 
and find some spherically symmetric exact solutions for /(i?)-gravity. In Sec IIV1 we review the Newman- Janis 
method to obtain axially symmetric solutions starting from spherically symmetric ones. The resulting metric is 
written in terms of two arbitrary functions. A further suitable coordinate transformation allows to write the metric 
in the so called Boyer-Lindquist coordinates. Such a transformation makes the physical interpretation much clearer 
and reduces the amount of algebra required to calculate the metric properties. In Sec|Vj the Newman- Janis method 
is applied to the spherically symmetric exact solution, previously derived by the Noether Symmetry, and an axially 
symmetric exact solution is obtained. This result shows that the Newman- Janis method works also in /(i?)-gravity. 
A physical application of the result is discussed in Sec IVIl Discussion ans concluding remarks are drawn in Sec. I VIII 

II. SPHERICAL SYMMETRY IN / (R) - GRAVITY 

Let us consider an analytic function f(R) of the Ricci scalar R in four dimensions. The variational principle for 
this action is: 



S / d 4 x s 



f(R) + XL r 



(1) 



where X — — j— , C m is the standard matter Lagrangian and g is the determinant of the metric 1 . 

c 4 

By varying with respect to the metric, we obtain the field equations 2 

= f'(R)R^ - \f{R)g^ - f(R).^ + g^nf(R) = XT^ 

(2) 

H = gP°U pa = 3Df'{R) + f(R)R-2f(R) = XT 



-2 J(V=g£ m ) 
~g 5g^ 

the trace. The most general spherically symmetric solution can be written as follows : 



where T pv — / f uii is the energy-momentum tensor of standard fluid matter and the second equation is 



We are adopting the convention Rp V = R p ' p, P v for the Ricci tensor and R a pp V = — for the Riemann tensor. Connections are 

Levi-Civita : 

V afj = ^9 tiP {9ap,fl + 9/3p,a ~ 9a?l,p) ■ 

2 It is possible to take into account also the Palatini approach in which the metric g and the connection T are considered independent 
variables (see for example [2lj1. Here we will consider the Levi-Civita connection and will use the metric approach. Sec [3, [22] for a 
detailed comparison between the two pictures. 
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ds z = m 1 (t',r')dt u +m 2 (t' ,r')dr u +m 3 (f ,r')dt'dr' + m A {t' ,r')dQ. 



(3) 



where m,, are functions of the radius r' and of the time t' . dfl is the solid angle. We can consider a coordinate 
transformation that maps the metric ([3]) in a new one where the off- diagonal term vanishes and m^t'^r') = —r 2 , 
that is 3 : 



ds 2 = gtt(t, r)dt 2 — g rr (t, r)dr 2 — r 2 dVl . 



(4) 



This expression can be considered, without loss of generality, as the most general definition of a spherically symmetric 
metric compatible with a pscudo - Riemannian manifold without torsion. Actually, by inserting this metric into the 
field Eqs. (O, one obtains: 



f'{R)R^ v — \ f{R)g l _ Ll , + H^u — XT^ 
f'(R)R-2f{R) +H = XT 



where the two quantities r H. ilv and H read 
R .... - 



'H-tiv — —f"(R)\ R,nv ~ r^-R t — r^„ij >r — g^v 



[g l \ t + g tt (In y/=g) ^ R tt + (g r \ r + g rr (In , 



R, 



+g u R.u + g rr R,rr 



-f"'(R) 



R,nR,v — g^v ( g tt R,t + g rr R, r 2 



n = g aT H aT = 3f"(R) 



^g tt >t + g u (in V=t)^R,t + U r , r + g rr (in V=g) >r ) R, r + g tt R <t t + g rr R, 



-3/'"(i?) 



g tt R/ + g rr R, r 2 



(5) 



(6) 



(7) 



Our task is now to find out exact spherically symmetric solutions. 

In the case of time-independent metric, i.e., gtt — a(r) and g rr = bjr), the Ricci scalar can be recast as a Bernoulli 
equation of index two with respect to the metric potential b(r) (see [9| for details) : 



h , {r) | , r 2 a\r) 2 - Aa(r) 2 -2ra{r)[2a{r)' +ra(r)"] \ , 



ra(r)[4a(r) +ra'(r)] 
where R = R(r) is the Ricci scalar. A general solution of © is: 



2 + r 2 R(r) 
4o(r) + ra'(r) 



b(r) 2 = 0. 



(8) 



b(r) 



exp[— J dr h(r)] 



K + fdrl(r) exp[- / dr h(r)] ' 



(9) 



where K is an integration constant while h(r) and l(r) are two functions that, according to Eq.(|8]), define the 
coefficients of the quadratic and the linear term with respect to b{r) [23|. We can fix l(r) = 0; this choice allows to 

2 

find out solutions with a Ricci scalar scaling as — - in term of the radial coordinate. On the other hand, it is not 

possible to have h(r) = since, otherwise, we get imaginary solutions. A particular consideration deserves the limit 
r — > oo. In order to achieve a gravitational potential b(r) with the correct Minkowski limit, both h(r) and l(r) have 
to go to zero at infinity, provided that the quantity r 2 R(r) turns out to be constant: this result implies b'(r) = 0, 
and, finally, also the metric potential b(r) has a correct Minkowski limit. 



This condition allows to obtain the standard definition of the circumference with radius r. 
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In general, if we ask for the asymptotic flatness of the metric as a feature of the theory, the Ricci scalar has to 
evolve to infinity as r~ n with n 2. Formally, it has to be: 

lim r 2 R(r) = r~ n , (10) 

r— ¥oo 

with n G N. Any other behavior of the Ricci scalar could affect the requirement to achieve a correct asymptotic 
flatness. 

The case of constant curvature is equivalent to GR with a cosmological constant and the solution is time independent. 



This result is well known (see, for example, |24() but we report, for the sake of completeness, some considerations 
related with it in order to deal with more general cases where a radial dependence for the Ricci scalar is supposed. If 
the scalar curvature is constant (R = Rq), field Eqs.([5]), being % tlv = 0, reduce to: 

(11) 

/oA) - 2/o = XT 

where f(Ro) = fo, /'(Rq) = fb- A general solution, when one considers a stress-energy tensor of perfect-fluid 
T MV = {p + p)u^u v -pg^v, is 



ds 2 = ( i + *i + ?*Rz±A a* dr \ - r Hn. (12) 

when p = —p. A = — -Mr and q^ 1 = f' Q . This result means that any /(i?)-model, in the case of constant curvature, 
exhibits solutions with de Sitter-like behavior. This is one of the reasons why the dark energy issue can be addressed 
using these theories [l| . 

If f(R) is analytic, it is possible to write the series: 

f(R) = A + * i? + *(i?), (13) 

where is a coupling constant, A plays the role of the cosmological constant and ^(i?) is a generic analytic function 
of R satisfying the condition 



lim R- 2 ^(R) = #1 , (14) 

where $i is a constant. If we neglect the cosmological constant A and is set to zero, we obtain a new class of 
theories which, in the limit R — > 0, does not reproduce GR (from Eq.tfM]), we have limi{_>o f{R) ~ R 2 )- I n such a 
case, analyzing the whole set of Eqs. (fTTj) . one can observe that both zero and constant ^ curvature solutions are 
possible. In particular, if R = Rq = field equations are solved for any form of gravitational potential entering the 
spherically symmetric background, provided that the Bernoulli Eq. ([5]), relating these functions, is fulfilled for the 
particular case R(r) = 0. The solutions are thus defined by the relation 



exp[— J dr h(r 



J r[a(r)-\-ra' (r)] 



being gtt(t, r) — b(r) from Eq. (J4j) . In some examples of /(-R)-models admitting solutions with constant^ or null 
scalar curvature are discussed. 



III. THE NOETHER SYMMETRY APPROACH AND THE SPHERICAL SYMMETRY 



Besides spherically symmetric solutions with constant curvature scalar, also solutions with the Ricci scalar depending 
on radial coordinate r are possible in /(i?)-gravity Q. Furthermore, spherically symmetric solutions can be achieved 
starting from a point-like /(i?)-Lagrangian [6(. Such a Lagrangian can be obtained by imposing the spherical symmetry 
directly in the action ([T]). As a consequence, the infinite number of degrees of freedom of the original field theory 
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will be reduced to a finite number. The technique is based on the choice of a suitable Lagrange multiplier defined by 
assuming the Ricci scalar, argument of the function f(R) in spherical symmetry. 

Starting from the above considerations, a static spherically symmetric metric can be expressed as 



ds 2 



and then the point-like f(R) Lagrangian 4 is 



A(r)dt 2 - B{r)dr 2 - M(r)dVL , 



(16) 



A 1 ' 2 f 9 f 
r = J —M' - 

2MB 1 / 2 A 1 / 2 B 1 / 2 



A'M 1 



Mf" 



A 1 / 2 B 1 / 2 



04 1/2 f" 

A 'R' — — R' M' - A 1 ' 2 B 1 ' 2 [(2 + MR)f - Mf] , (17) 



which is canonical since only the configuration variables and their first order derivatives with respect to the radial 
coordinate r are present. Details of calculations are in Q. Eq. (|17j) can be recast in a more compact form introducing 
the matrix representation : 



C = q ' fq' + V 



(18) 



where q = (A,B,M,R) and q' = (A' , B' , M' , R') are the generalized positions and velocities associated to C. It is 
easy to check the complete analogy between the field equation approach and point- like Lagrangian approach [|| . 

In order to find out solutions for the Lagrangian (JT7J), we can search for symmetries related to cyclic variables and 
then reduce dynamics. This approach allows, in principle, to select /(i?) -gravity models compatible with spherical 
symmetry. As a general remark, the Noether Theorem states that conserved quantities are related to the existence of 
cyclic variables into dynamics [25l - l27| . 

d 2 C 

It is worth noticing that the Hessian determinant of Eq. (|17l) , 



dq'idq'j 



, is zero. This result clearly depends on the 



absence of the generalized velocity B' into the point -like Lagrangian. As matter of fact, using a point-like Lagrangian 
approach implies that the metric variable B does not contributes to dynamics, but the equation of motion for B has 
to be considered as a further constraint equation. Then the Lagrangian (|17p has three degrees of freedom and not 
four, as one should expect a priori. 

Now, since the equation of motion describing the evolution of the metric potential B does not depend on its 
derivative, it can be explicitly solved in term of B as a function of the other coordinates : 



B = 



2M 2 f"A'R' + IMf'A'M' + AAMf'M'R' + Af'M'' 
2AM[(2 + MR)f' - Mf] 



(19) 



By inserting Eq. (|19j) into the Lagrangian (|17j). we obtain a non- vanishing Hessian matrix removing the singular 
dynamics. The new Lagrangian reads 5 



with 



L = q' tiq' = 



[(2 + MR)f - f M] , 
M 



[2M 2 f"A'R' + 2MM\f'A' + 2Af"R!) + Af'M' 2 ] 



(20) 



If one assumes the spherical symmetry, the role of the afftne parameter is played by the coordinate radius r. In this 
case, the configuration space is given by Q = {A, M, R} and the tangent space by TQ = {A, A' , M, M', R, R'}. On 
the other hand, according to the Noether Theorem, the existence of a symmetry for dynamics described by Lagrangian 
(fT7| implies a constant of motion. Let us apply the Lie derivative to the (fT7|) , we have 6 : 



L X L = a ■ V„L + ol ■ V„'L 



a ■ V 9 L + 2 V 9 a 



(21) 



4 Obviously, the above choices are recovered for A(r) = a(r), B{r) = b(r), and M(r) = r 2 . Here we deal with A,B,M as a set of 
coordinates in a configuration space. 

5 Lowering the dimension of configuration space through the substitution I I19H does not affect the dynamics since B is a non-evolving 
quantity. In fact, inserting Eq. I|19|l into the dynamical equations given by J2J, they coincide with those derived by (JTtJ . 

6 From now on, q indicates the vector {A, M, R}. 
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that vanishes if the functions a satisfy the following system 

a ■ V 9 L + 2(V 9 a)*L = — > a, + 2^L 4m = . (22) 

Solving the system (f2"2"j) means to find out the functions which assign the Noether vector [25|, [2(| . However the 
system (|22[) implicitly depends on the form of /(i?) and then, by solving it, we get also /(i?)-models compatible with 
spherical symmetry. On the other hand, by choosing the /(i?)-form, we can explicitly solve (|22[) . As an example, one 
finds that the system (|2"2"|) is satisfied if we choose 



f(R) = f R s , a= (ai,a 2 ,a 3 ) = ( (3 - 2s)kA, -kM, kRj , (23) 

with s a real number, k an integration constant and /q a dimensional coupling constant 7 . This means that, for any 
f(R) = R s , exists, at least, a Noether symmetry and a related constant of motion S : 

So = a ■ Vq'L = 2skMR 2s ^[2s + (s - l)MR][(s - 2)RA' - (2s 2 - 3s + 1)AR'] . 

A physical interpretation of So is possible if one gives an interpretation of this quantity in GR, that means for 
f(R) = R and s = 1. In other words, the above procedure has to be applied to the Lagrangian of GR. We obtain the 
solution 

a GR = (-kA, kM) . (24) 
The functions A and M give the Schwarzschild solution, and then the constant of motion acquires the standard form 

v 2GM ton 
S = — J - . (25) 

<r 

In other words, in the case of Einstein gravity, the Noether symmetry gives, as a conserved quantity, the Schwarzschild 
radius or the mass of the gravitating system. This result can be assumed as a consistency check. 
In the general case, f(R) = R s , the Lagrangian (fl7|) becomes 

sR 2s - 3 \2s + (s-l)MR\ r , , / , ,,,,■>, 

L = ■ — --j -[2{s - l)M 2 A'R' + 2MRM A + 4(s - I) AM MR + ARM' 2 } , 

and the expression (fTTj)) for B is 

_ s[2(s - l)M 2 A'R' + 2MRM'A' + 4(s - \)AMM'R! + ARM' 2 } 

~ 2AMR[2s + (s - 1)MR] ^ ' 

As it can be easily checked, GR is recovered for s = 1. 

Using the constant of motion (j2"4"]l . we solve in term of A and obtain 



2s A -35 + 1 



4s* -9s + 5 



A = R ~ 1 kl + S f ; R p— — ^ r A (27) 

\ U J 2ks(s - 2)M[2s + (s - l)MR} J v ; 
for s ^ 2 and fci an integration constant. For s = 2, one finds 

4 ^2 (OR) 

12kr 2 (4 + r 2 R)RR' ' 1 ; 

These relations allow to find out general solutions for the field equations giving the dependence of the Ricci scalar on 
the radial coordinate r. For example, a solution is found for 

s = 5/4 , M = r 2 , R = 5r~ 2 , (29) 



7 The dimensions are given by R 1 s in terms of the Ricci scalar. For the sake of simplicity, we will put /rj = 1 in the forthcoming 
discussion. 



7 



obtaining the spherically symmetric space-time 

ds 2 = (a + (3r)dt 2 - l „ dr 2 - r 2 dQ , (30) 
2 a + pr 

where a is a combination of T, and k and /3 = ki. In principle, the same procedure can be worked out any time Noether 
symmetries are identified. Our task is now to show how, from a spherically symmetric solution, one can generate an 
axially symmetric solution adopting the Newman- Janis procedure that works in GR. In general, the approach is not 
immediately straightforward since, as soon as f(R) 7^ R, we are dealing with fourth-order field equations which have, 
in principle, different existence theorems and boundary conditions. However, the existence of the Noether symmetry 
guarantees the consistency of the chosen /(i?)-model with the field equations. 



IV. AXIAL SYMMETRY DERIVED FROM SPHERICALLY SYMMETRIC SOLUTIONS 



We want to show now how it is possible to obtain an axially symmetric solution starting from a spherically symmetric 
one adopting the method developed by Newman and Janis in GR. Such an algorithm can be applied to a static 
spherically symmetric metric considered as a "seed" metric. Let us recast the spherically symmetric metric (|4]) in the 
form 



ds 2 = e 2 ^dt 2 - e 2X ^dr 2 - r 2 dn, 



(31) 



with g tt (t,r) = e 2 ^ and g rr (t,r) = e 2X{ - r \ Such a for m is suitable for the considerations below. Following Newman 
and Janis, Eq. (|3"Tj) can be written in the so called Eddington-Finkelstein coordinates (u, r, 0, </>), i.e. the g rr component 
is eliminated by a change of coordinates and a cross term is introduced |28j . Specifically this is achieved by defining 
the time coordinate as dt — du + F(r)dr and setting F(r) = ±e x ^~^^ r \ Once such a transformation is performed, 
the metric (f3"Tj) becomes 



ds 2 = 



^du 2 ±2e x ^ + ^dudr 



r 2 dn. 



(32) 



The surface u = costant is a light cone starting from the origin r = 0. The metric tensor for the line element (|32p 
in null-coordinates is 



/ 

± e -A(r) 







-4>{r) 



_ e -2A(r) 






The matrix (|33[) can be written in terms of a null tetrad as 






-1/r 2 








-l/(r 2 sin 2 6>) / 



(33) 



a 



Vn" + rn" - rn^fh" - m"fh ti , 



(34) 



where n^, and m M are the vectors satisfying the conditions 



I pi 1 * = m^in^ = n^n^ = 0, 1^ = -m^m^ = 1, l^m^ = n fJ ,m fl = 0. (35) 

The bar indicates the complex conjugation. At any point in space, the tetrad can be chosen in the following manner: 

is the outward null vector tangent to the cone, n M is the inward null vector pointing toward the origin, and 
and m' 1 are the vectors tangent to the two-dimensional sphere defined by constant r and u. For the spacetime (f33j) . 
the tetrad null vectors can be 



(36) 
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Now wc need to extend the set of coordinates x^ — (u, r, 9, (f) replacing the real radial coordinate by a complex 
variable. Then the tetrad null vectors become 8 



A new metric is obtained by making a complex coordinates transformation 



(37) 



a? -> - +ij/ M (a: <T ). 



(38) 



where y M (x cr ) are analityc functions of the real coordinates x a , and simultaneously let the null tetrad vectors Z£ 
(Z M , , m/ 1 , m M ) , with a = 1,2,3,4, undergo the transformation 



(39) 



Obviously, one has to recover the old tetrads and metric as soon as x a = x a . In summary, the effect of the " tilde 
transformation" (|38l) is to generate a new metric whose components are (real) functions of complex variables, that is 



with 



For our aims, we can make the choice 



9nu ^ dins -xxx i y 



' u = u + ia cos ( 



x" = x" + ia{6f -6%)cosl 



r = r — ia cos t 



= e 



where a is constant and the tetrad null vectors (|37[) . if we choose r = r, become 



(40) 



(41) 



(42) 



l" = 5% 



= _i e -2A(f,e)^ ■ + e -\(f,o)-<K?,o) s * 



rnr 



*/2{r—ia cos 9) 



\^2(r+ia cos 9) 



ia(<yff-^) sin (9 + ^ + ^ 



(43) 



It is worth noticing that a certain arbitrariness is present in the complcxification process of the functions A and </>. Obviously, we have 
to obtain the metric J33D as soon as r = f . 
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From the transformed null tetrad vectors, a new metric is recovered using (|34|) . For the null tetrad vectors given 
by (|43[) and the transformation given by ((32), the new metric, with coordinates x' 1 = (u, f, 9, 0), is 



/ a 2 sin 2 6 „-A(f ,0)-6(r ,6) , a 2 sin 2 
/ y2 e ~r ^2 



_ p — 2\(r : 6) a~ sin 



V 



o 
o 

i 

S 2 



--2- \ 
£ 2 \ 







(44) 



i / 

S 2 sin 2 fl / 



where £ = V^ 2 + a 2 cos 2 6*. In the covariant form, the metric (|4"4"() is 



/ e 20(f=,e) e A(f,e)+0(f,e) 




V 



o 



-s 2 



-[£ 2 



a 2 sin 2 



a( >(^)+A(r,0) sill 2 


A(f,6) _ 0(f,0)M„:„2 



>(2e> 



(45) 



)] sin 2 9 ) 



Since the metric is symmetric, the dots in the matrix are used to indicate = g u ^. The form of this metric gives 
the general result of the Newman- Janis algorithm starting from any spherically symmetric "seed" metric. 

The metric given in Eq. (|4"5"j) can be simplified by a further gauge transformation so that the only off-diagonal 
component is g^f This procedure makes it easier to compare with the standard Boycr-Lindquist form of the Kerr 
metric [28| and to interpret physical properties such as the frame dragging. The coordinates u and tfi can be redefined 
in such a way that the metric in the new coordinate system has the properties described above. More explicitly, if we 
define the coordinates in the following way 



du = dt + g(f)dr and d<p — d(j) + h(f)dr 



(46) 



where 



9(f) 



S A ^- e >(S 2 +a 2 sin 2 8e^'<» + <f^'<») 

e *(f,8)( S 2 +a 2 sin 2 ee 2A(f,8)) 



I V I — S 2 +a 2 sin 2 9e 2A(f,8) 

after some algebraic manipulations, one finds that, in (t,r,0,<j>) coordinates system, the metric (|45[) becomes 

\ 



(47) 



9^ 



( e 2 ^ 

-E 2 /(S 2 e- 2A ^) +a 2 sin 2 0) 



V 



-p 2 






a 2 sin 2 e e ^e)( 2e Hm „ e <t>{r,e))] sin 2 9 j 



(48) 



This metric represents the complete family of metrics that may be obtained by performing the Newman- Janis 
algorithm on any static spherically symmetric "seed" metric, written in Boycr-Lindquist type coordinates. The 
validity of these transformations requires the condition S 2 + a 2 sin 

e 2A(f,0) o, where e 2Xl ~ f ^ > 0. Our task is now 
to show that such an approach can be used to derive axially symmetric solutions also in /(i?)-gravity. 

V. AXIALLY SYMMETRIC SOLUTIONS IN /(TJ)-GRAVITY: AN EXAMPLE 



Starting from the above spherically symmetric solution (|30p . the metric tensor, written in the Eddington-Finkelstein 
coordinates (u, r, 0, </>) of the form ((551) is 



(0 



\ ■ 



2_ 

fir 

_n _ 2a 

fir 






-1/r 2 

-l/(r 2 sin 2 i 



(49) 
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The complex tetrad null vectors (|37[) are now 



°1 + V /3 W 



y^f V 2 sin J / 

By computing the complex coordinates transformation (|42[) . the tetrad null vectors become 



(50) 



q Rc{f} 



/9 



Vs°o 



(51) 



Now by performing the same procedure as in previous section, we derive an axially symmetric metric of the form 
but starting from the spherically symmetric metric (|30[) . that is 



9nv 



I r(a+f3r)+a 2 cos 2 6 



V 



o 

/3S 2 






-E 2 



q(-2ar-2 ) 3S 2 +V2?S 3/2 ) sin 2 6 
22 




^2 a 2 (ar+ffS 2 -V2?S 3/2 )sin 2 6 
J 2 



sin 



\ 



/ 



(52) 



It is worth noticing that the condition a = immediately gives the metric (|30[) . This is nothing else but an example: 
the method is general and can be extended to any spherically symmetric solution derived in /(i?)-gravity. 



VI. A PHYSICAL APPLICATION: GEODESICS AND ORBITS 

Let us discuss now a physical application of the above result. We will take into account a freely falling particle 
moving in the space-time described by the metric (|52j) . For our aims, we make explicit use of the Hamiltonian 
formalism. Given a metric g^, the motion along the geodesies can is described by the Lagrangian 



C{x^x») = -g^x" 



(53) 



where the ovcrdot stands for derivative with respect to an affine parameter A used to parametrize the curve. The 
Hamiltonian description is achieved by considering the canonical momenta and the Hamiltonian function 



dC 

dx^ 



n = - c , 



(54) 



1 



that results TL = —PftPv9^ u - The advantage of the Hamiltonian formalism with respect to the Lagrangian one is that 
the resulting equations of motion do not contain any sign ambiguity coming from turning points in the orbits (see, 
for example, j29j ) . The Hamiltonian results explicitly independent of time and it is 7i = ——m 2 , where the rest mass 
to is a constant (to = for photons). The geodesic equations are 

dx» dH 



dX <9p M 



q^p^p", 



(55) 



dX 



&H_ 

dxfi 



a/3 



I dg 

2 dx» 



■PaPfi = g^FZ-yPaPp 



(56) 
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In addition, since the Hamiltonian is independent of the afnnc parameter A, one can directly use the coordinate time 
as integration parameter. The problem is so reduced to solve six equations of motion. Using the above definitions, it 
is easy to achieve the reduced Hamiltonian (now linear in the momenta) 



H = -pa = 



Pi9 



Oi 



.,00 



Pi 9 



Qi 



-.00 



m 2 + PiPjg 13 
„oo 



1/2 



(57) 



with the equations of motion 



dx l _ dH 
dt dpi ' 



dpi dH 
dt dx l 



(58) 



that give the orbits. The method can be applied to the above solution (|52[) considering the following line element 



ds- 



r{a + j3r) + a 2 /3 cos 2 2 a{-2ar - 2f3Y? + V2/3E 3 / 2 ) sin 2 



E 

/3E 2 



-dt 2 + 2 



-dr 2 - Yrd0 



2E 



dtd<fr 



l {ar + /3E 2 - v^E 3 / 2 ) sin 2 



2ar + ^(a 2 +r 2 + E 2 ) 
by which the elements of the inverse metric can be easily obtained: 



sm 



2 fl ji2 



4E 2 



2 a 2 sin 2 9(rQ-v / 2v / ?S 3/2 +/3S 2 ) 

— ^ 



i 2 sin (2ra - v^v^S 3 / 2 + 2£E 2 ) + 4E (a 2 f3 cos 2 + r(r/3 + a)) ( E 2 
(3 (a 2 + r 2 + E 2 ) + 2ra 



a 2 sin 2 9(ra-v / 2 v / ?S 3 / 2 + / 3S 2 ) 



/?E 2 



1 

E 2 



.9 = -^r 



, ttj> 



2aE (-2m + ^V?E 3 / 2 - 2/3E 2 



-■ 2 sin 2 (2ra - v^v^S 3 / 2 + 2^E 2 ) Z + 4E [a 2 /3 cos 2 + r(r/3 + a)} 



2 fl2 sin 2 # (Va — 

E ?i 



4E esc 2 [a 2 /3 cos 2 + r(r/3 + a)] 



a 2 sin 2 (2ra - \/2\/?E 3 / 2 + 2/3E 2 ) 2 + 4E [a 2 (3 cos 2 + r(r/3 + a)} 



E 2 



a 2 sin 2 e(ra-V2 v / FE 3 / 2 +,9E 2 ) 



and the null ones 



tr td rd rd> 6cb n 

9 =9 =9 =9~9=0. 



(59) 



(60) 



Let us consider the equatorial plane, i.e. 0=^,0 = 0, and assume a = 1 and /3 = 2. The reduced Hamiltonian 
can be written as 



H{r,0,(t>,p r ,Pe-,Pcj>;t) 



2 a p, (-2r3 + r=» -l) r r , _ s 

a 2 (-2(r - l)r 2 - 1) + r 5 U\ 1 ? K ' 



-a 2 (-2(r - l)r 2 - l) - r 5 ) (a 2 (r 2 (r{2r - 3)(2r + 1) + 6) - 2) + (2r + l)r 4 ) x 
lty(2r+l) Pr (a 2 +r 2 +r) + Pe 



a 2 (r 2 {r(2r - 3)(2r + 1) + 6) - 2) + (2r + l)r 4 



-p r + l 



(61) 
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270 



Figure 1: Relative motion of the test particle with m = 1. 

It is independent of <j> (i.e. we are considering an azimuthally symmetric spacetime), and then the conjugate momentum 
is an integral of motion. From Eqs. (|58p. one can derive the coupled equations for {r,0,<fi,p r , pg} and integrate 
them numerically (the expressions are very cumbersome and will not be reported here). To this goal, we have to 
specify the initial value of the position-momentum vector in the phase space. A Runge-Kutta method can be used to 
solve the differential equations. In Fig[TJ the relative trajectories are sketched. 

VII. DISCUSSION AND CONCLUDING REMARKS 

We have shown that the Newman- Janis method, used to derive axially symmetric solutions in GR, works also in 
/(-R)-gravity. In principle, it could be consistently applied any time a spherically symmetric solutions is derived. The 
method docs not depend on the field equations but directly works on the solutions that, a posteriori, has to be checked 
to fulfill the field equations. 

The key point of the method is to find out a suitable complex transformation which, from a physical viewpoint, 
corresponds to the fact that we are reducing the number of independent Killing vectors. From a mathematical 
viewpoint, it is useful since allows to overcome the problem of a direct search for axially symmetric solutions that, 
in /(i?)-gravity, could be extremely cumbersome due to the fourth-order field equations. However, other generating 
techniques exist and all of them should be explored in order to completely extend solutions of GR to /(ii)-gravity. 
They can be more general and solid than the Newman- Janis approach. A good source for references and basic 
features of generating techniques is reference (T^j. In particular, the paper by Talbot [l4|, considering the Newman- 
Penrose approach to twisting degenerate metrics, provides some theoretical justification for the scope and limitations 
of adopting the "complex trick". As reported in Chap. 21 of [l3j |. several techniques can be pursued to achieve 
axially symmetric solutions which can be particularly useful to deal with non-empty space-times ( in particular when 
perfect fluids are the sources of the field equations) and to deal, in general, with problems related to Einstein-Maxwell 
field equations. We have to stress that the utility of generating techniques is not simply to obtain a new metric, 
but a metric of a new spacetime with specific properties as the transformation properties of the energy-momentum 
tensor and Killing vectors. In its original application, the Newman- Janis procedure transforms an Einstein-Maxwell 
solution (Reissner- Nordstrom) into another Einstein-Maxwell solution (Kerr- Newman). As a particular case (setting 
the charge to zero) it is possible to achieve the transformation between two vacuum solutions (Schwarzschild and 
Kerr). Also in case of /(i?)-gravity, new features emerge by adopting such a technique. In particular, it is worth 
studying how certain features of spherically simmetric metrics, derived in /(i?)-gravity, result transformed in the 
new axially symmetric solutions. For example, considering the f(R) spherically symmetric solution studied here, the 
Ricci scalar evolves as r~ 2 and then the asymptotic flatness is recovered. Let us consider now the axially symmetric 
metric achieved by the Newman- Janis method. The parameter a ^ indicates that the spherical symmetry (a = 0) 
is broken. Such a parameter can be immediately related to the presence of an axis of symmetry and then to the fact 
that a Killing vector, related to the angle 9, has been lost. To conclude, we can say that once the vacuum case is 
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discussed, more general spherical metrics can be transformed in new axially symmetric metrics adopting more general 
techniques [ToT ]. These approaches will be examined and discussed in future works. 
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